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This paper is devoted to a prediction of new effect - the occurrence of a heat flow perpendicular 
both to temperature gradient and magnetic field in molecular dielectric where rotary degrees of 
freedom of molecules are defrozen. The method of the moments considering processes of phonon 
scattering on molecules with their rotary condition changing is developed. 
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^ ; l.INTRODUCTION. 

Q-i. Recently there was a message [H, [3] about detection at presence of a magnetic field in ionic dielectric Tb 3 Ga 5 0i 2 a 
' heat flow perpendicular both to temperature gradient and magnetic field. This phenomenon is related to Righi-Leduc 
effect in metals where it is also observed a component of a heat flow, directed on [B, VT]. As in dielectrics there are no 
free charge carriers the physical reason of the phenomenon is far from Hall effect and Righi-Leduc effect mechanisms. 
In works [H], (3] and theoretical works [3] the observable effect communicates from spin-orbital phonon interaction 
with magnetic biased spins of paramagnetic ions. However, there is another mechanism which can lead to occurrence 
_C! of transverse heat flow. In molecular gases this is an anisotropic scattering of molecules which rotary moments precess 
in external magnetic field B. The observation of a heat flow perpendicular both to temperature gradient and magnetic 
field in molecular gases is called odd Senftleben-Beenakker effect (jj. In the present work this effect is generalized 
on a case of molecular crystals where oscillating near lattice points molecules can rotate freely at the same time if 
temperature is above the rotary degrees of freedom freezing temperature [f| . Thus as translational motion of molecules 
in a solid body is absent, there is no direct analogue of classical Senftleben-Beenakker effect. Phonons transfer energy 
in dielectrics. In a considered solid body we assume there are no paramagnetic particles and spin-orbital mechanism 
of a magnetic field influence on heat conductivity is absent. However, if the molecule rotary moments of this body 
can rotate freely in a magnetic field these moments precess. Precession affects the probability of anisotropic phonon 
scattering at nonspherical molecules. It should lead to phonon energy flow change appeared in presence of temperature 
gradient and to occurrence of magnetic field dependence of thermal conductivity tensor: 
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00 ! H ik {B) = K ki {-B). (1) 

In the present work it is predicted, that the thermal Hall effect can be observed not only in ionic, but also in 
molecular dielectrics in the field of temperatures where rotary degrees of molecules are defrozen. Rotary moments 
M will be considered as classical vectors. (More realistic case of quantum rotation will be investigated separately). 
The considered phenomenon is caused by correlation of phonon movement and the rotary moments and described by 
joint distribution function of phonon and rotary degrees of freedom / (p, Mj . Correlation exists, if phonon scattering 

at the rotary moments is anisotropic. For an estimation of such anisotropy it is necessary to solve a mechanical 
problem about interaction of two oscillating near adjacent lattice points molecules. We shall not do this, and using 
a method of the moments for solution of the corresponding kinetic equation the eigenvalues of the collision operator 
we shall consider as adjustable parameters. Thus, three parameters are entered into theories instead of one usual 
adjustable parameter (phonon relaxation time r): 1) reciprocal phonon lifetime owing to anharmonicity or dispersion 
on defects Qq = 1/r, 2) phonon collision frequency with change of a rotary condition of molecules f2i = AS1 (A< 1 
- nonsphericalness parameter) and 3) precession frequency lub = — 7-B. (7 - gyromagnetic ratio). It is used spherical 
coordinates and system of units with H = 1, ks = 1 here. It will be shown, that the antisymmetric part of heat 
conduction tensor is proportional to a magnetic field 

xT(B)~\ 2 (u B /n). (2) 

2. THE KINETIC EQUATION. 

So, we shall find a magnetic field influence on heat conduction of a molecular crystal with rotary degrees of freedom. 
We investigate a simple model of substance in which all molecules or their part have the rotary moments. We shall 
consider system as a solution of Debye phonon gas with a spectrum ui p — cp, and located in units of lattice of two-atom 
molecules classical rotary moments M noninerating with each other. (In contrast to a problem of heat conduction in 
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a paramagnetic with spin-phonon interaction Q phonon polarization does not play a role and detailed elaboration 
of acoustic phonon properties does not matter). Phonon and rotary moments interaction is not enough, and in 
gas approximation we shall describe evolution of system by Boltzmann equation [?J for joint distribution function 
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+ Stf = 0. (3) 



Here f , v = du/dp are a phonon coordinate and velocity (\v\ = c). If the subsystem of rotary moments could move like 
gas in space, it would be necessary to bring in the left part of J3]) an expression ^^-gp-> where fu is a coordinate of a 
point where one of the rotary moments is located. But such motion is absent and = 0. The fourth term at the left 
describes rotary moment precession in a magnetic field (jfiM = 7 [Mi?] J . In paramagnetic substance, there is a mix 
of components with different electron spin directions. The odd effect is created by component, consisting of molecules 
which have a spin projection on molecule axis equal to zero. In nonparamagnetic body molecule gyromagnetic ratio 7 
is much less, but not zero because of nuclei slipping in relation to electron shell [§]. For simplicity we shall consider a 
case of nonparamagnetic molecules. Below we shall generally use spherical coordinate representation where precession 
operator has as much simple form as possible: 

(l B = j[MB}— = uj B - , uj b = -7-B, n B Y hm2 (M) = im 2 L0 B Yi 2m2 (M) (4) 

dM d( PM 

In a state of thermodynamic equilibrium we have 

/W(V,M) =N^(p)^°\M), (5) 

where 

7V(°)(p) = (ex P K/T)-l)- 1 , ^)(M) = -J_exp(-e/T) (6) 

(e = M 2 /2I and I are energy and the moment of inertia of a two-atom molecule). The phonon number is set by 
temperature, and the rotary moment number is fixed. Distribution function (tp^(M)) is referred to one molecule. 
At the decision of a stationary problem of heat conductivity the first term in |(3]) is absent, and the second is equal to 

(tfV)/<°> = (vVN^(p,r)) ip<® = (vVT) {uo/T 2 )N^ (n<® + l) ^ = / (0) QVT, Q = v^((n^ + l) (7) 

in linear approximation on a temperature gradient. The rotary moments are attached to a lattice and cannot give 
the contribution to a heat flow. Therefore calculating a gradient we have neglected spatial heterogeneity. 
Let's represent distribution function in the form of 

/ = /(°)(l+xVT) (8) 

In linear approximation on VT the collision integral is expressed through the linear integral operator 

St/ = (VT)/(°)fix, (9) 

which acts on x- Positive definite Hermitian collision operator Cl is completely described by the matrix elements 
Qab = (ip*£lipbj- Matrix elements fl a b possess properties: 

(iPZWb)* = (VWa) = ((fi^a) r b ) , (r a Wa) > (10) 

As it has been already noted in introduction, we shall not write explicit rotary moment dependence of phonon 
scattering cross-section and we shall apply standard approach in kinetics when dispersion properties are described 
phenomenologically by il a b which are assumed known, i.e., actually, they are adjustable parameters. Thus, removing 
from J3]) a time derivative, using (iTj) , , (J9j , and reducing on /'°'VT, we obtain the vector integro-differential equation. 

(n + n B )x + Q = o, (11) 
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3. METHOD OF MOMENTS. 

Value of a heat flow and heat conduction tensor in Debye model are defined by 

Ql = 5> iW )/ = / ^3 M (/ (0) xVr) = -K ik (VT) k , K lk = - ( Vl oj Xk ) , (12) 

where (A) = J ^ d -^f(o) A . 

So, the problem comes to the solving of vector equation l[TTj) . Vector-function x we shall search in the form of 
expansion 
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YsXia^a^M) (13) 



on orthonormal vector functions ip a (p, M^j : (ip*ipb) — $ab- Orthonormal functions ip a are convenient for choosing 
in the form of superposition of products of spherical polynoms depending on impulse directions Y/ imi (p) and rotary 
moment Yi 2m2 (M) 

^Wxfa ~ E C M 2 m 2 Y kmi (p)Y h m 2 (M), (14) 



mi m,2 



where C; 1 1 ™ 1 ; 2 „ l2 is Clebsch-Gordan coefficients. To simplify notations, we do not write explicitly normalization 
and factors containing scalar polynoms depending on translational and rotary energies. Transition from cartesian 

coordinates to spherical is carried out by coefficients ti m , where t xm = , t ym = 1 -j\J~^§-, m = ±1 

Qi = E ^-mQlm, *ik = ~ (viUJXk) = - E t *im t km {v* lm UJXl m ) (15) 
m tci 

In particular, off-diagonal element of heat conduction tensor x yx is proportional to Im (t) u wxn). It is responsible for 
a heat flow in a direction (y) , perpendicular to temperature gradient (x) and magnetic field (z) . For solving equation 
(fTTTl in 2-moment approximation as a first function from lfT3j) we shall choose wave function contained in Q , which 
describes heterogeneity in l[TTj) . i.e. 

iplmio = v lmU (16) 

As a second function we shall choose vector-function which is even on M and odd on p. The simplest expression, 
satisfying to this requirement in cartesian coordinates is M or in spherical coordinates 

VW2= E ^Im^^lmx^WM)- (17) 
mim2 

Term Im(w u wxu) is easy to find, if the probability of molecula rotary condition changing is small (weak non- 
sphericity) at phonon dispersion. 

Let in zeroth-order approximation the collision operator in representation of functions lfl4|l be reduced by one 
relaxation time 

ft< >V„ = no^», no = i/T, (is) 

and off-diagonal part AfiW, proportional to small parameter (A <C 1), will transform the functions depending only 
from an impulse, to functions of a kind lfl4|) . depending from M 

O (1) ^i m io = CliipimU (19) 

The equation ([11]) 

(OW +Af2 (1) +n B ) X n +Qn -0 (20) 
is solved by decomposition of function x on degrees of A. We shall determine operator K by expression 

K = + he)- 1 (21) 



The equation ([20)1 we shall write down in the form of 

{k- 1 + xn^)xn + Qn = (22) 

Then we shall increase it on K at the left and then - on (1 + XKfl^)^ 1 . From here we find 

xii = -(i + xkCi^^kQh = -(i - xk&w + x 2 kn^kn^ + ...)kq u (23) 

The required term is equal to 

Im (v^uxu) = - Im (^^(1 - XKCl {1) + X 2 kCl w kCl {1) + ...)KQ n ) (24) 



As function Qn does not depend from M operator £Ib does not act on it, and 

KQ n = (OW)- x Qii = Qn/Qo. (25) 

As a result the right side 11241) becomes 

- Im /vj x w(l - A-^-fiW + X 2 -^Q^KQ^ + ...)Q n ) -j- (26) 

\ "0 "0 / "0 

First two terms give zero. The nonvanishing contribution gives a term of the second order on A, and accurate within 
numerical factor the off-diagonal part of heat conduction tensor is equal to 



^yx r>2 
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ClWMW ( VllUJ ) (n^ + 1^ (27) 



As the given frequencies integral quickly converges both at low and at high frequencies, then it is possible to neglect 
(iV(°) + l) and write in a compact kind 

Kyx ~^ Im ( P * kP )> (28) 
where function P = fiW (vhuj) is introduced. In concordance with l|16|l and (fl9| it is written in form of 

P = fii</W2- (29) 
To calculate a matrix element (^P*KP\ it is necessary to find function 

F = KP=(ClW +(l B )- 1 P (30) 
For this purpose we should solve following equation 

(m + «>b-£-) F = !li^ C]} ni2m2 Y lmi (p)Y 2m2 (M), (31) 

The solution will be searched in the form of 

F= J2 Fm 2 CH l2m2 Y lmi (f)Y 2m2 (M). (32) 

mi?7i2 

Expansion coefficients are easy to be found in approximation of one relaxation time: Cl^F — QqF. Then 

Ui +u> b -^-\f= (^o + im 2 u; B ) F ni2 C^ l2m2 Y lmi (p)Y 2m2 (M) = ^ £ C\^ ll2m2 Y\ mi (p)Y 2m2 (M) (33) 

^ r / mi m2 mim2 

From here we obtain 

F ni2 = fii (fi + im^gy 1 , (34) 
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and the value interesting us (flOjl becomes 
A 2 



Im 1 



/ V V 



E C lU2mi Y lm{ ® Y 2mi ( Al ) E ^ (°° + Wfl)" 1 C l2m /l mi ^ (M) ) (35) 



After averaging on p and M we come to final expression: 

A 2 Sli x - , n N 2 m 2 uJB 



mim2 



E ( c 

(36) 



i2 o m , ma ng + (m 2 w B ) 



We see, that the perpendicular heat flow is odd on a field and is maximal, when order of magnitude of £ = wbM 
is 1. 

Here the case of nonparamagnetic molecules was considered. In a solid body at achievable fields effect seems to be 
much less than the maximal value, and the perpendicular heat flow grows linearly with a field 

>c yx ~ A 2 o;b^o 4 (37) 
The even effect in a magnetic field is described by the valid part of expression which imaginary part was written in 

Xxx ~ —^r ^2 ( C lLi2m 2 ) o2 ; , 72 ( 38 ) 

U ° ™ n 2 Q + {m 2 LJ B ) 

Even and odd effects for gas were experimentally observed many times (see for example [H] and For a solid 
body with rotary degrees of freedom the of magnetic field influence on heat conductivity was not experimentally 
researched yet. 

In the present work it is shown, that in molecular crystals where quasi-free molecule rotation exists in wide range 
of temperatures, the heat conductivity effect, caused by occurrence of a heat flow perpendicular to a temperature 
gradient, similar to the Hall effect in metals should be observed odd on a magnetic field. 
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